Trapped ion emulation of electric dipole moment of neutral relativistic particles 
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The electric dipole moments of various neutral elementary particles, such as neutron, neutrinos, 
certain hypothetical dark matter particles and others, are predicted to exist by the standard model 
of high energy physics and various extensions of it. However, the predicted values are beyond the 
present experimental capabilities. We propose to simulate and emulate the electric dipole moment of 
neutral relativistic particles and the ensuing effects in the presence of electrostatic field by emulation 
of an extended Dirac equation in ion traps. 

PACS numbers: 11.30.Er 



The standard model of particle physics as well as all 
SUSY models predict that elementary particles should 
possess an intrinsic electric dipole moment (EDM) 
aligned with the spin of the particle. There are certain 
models where dark matter is thought to consist of electri- 
cally neutral relativistic particles possessing EDM and/or 
magnetic dipole moment (MDM) [l[. The existence of 
EDM generally requires the breaking of parity and time 
reversal symmetries. The validity of the CPT theorem 
then implies the breaking of CP symmetry. Mechanisms 
violating CP symmetry in the standard model are respon- 
sible for the generation of EDM of elementary particles 
with predicted values d n ~ 10 _32 e cm for the neutron 
and d e < 10~ 33 e cm for the electron 0. 

The CP violation is also connected with the problem 
of the observed matter-antimatter asymmetry in the uni- 
verse. The SUSY models, which address some of the 
shortcomings of the standard model, predict much larger 



values for the EDM of elementary particles: d, 

-25 



< 10" 



cm and 10~ 28 < d n < lCT 25 e cm Q. 

The experimental search for the EDM of the neutron, 
whose neutral charge makes it most suitable for EDM 
measurement, has been initiated some 60 years ago by 
Smith, Pureed and Ramsey Q. Decades of effort have 
also been dedicated to the measurement of the electron 
EDM. A long line of experiments spanning over 60 years 
have set only upper limits for the EDM's [1, \d e \ < 
1.05 x 10~ 27 e cm and \d n \ < 2.9 x 10~ 26 e cm. 

Recently, the field of quantum simulations of physical 
systems in ion traps has gained momentum including the 
simulation of relativistic systems. Proposals have been 
put forward for simulation of the free Dirac equation [?| , 
the Dirac oscillator || , the ID Dirac equation with vari- 
ous Poincare invariant potentials relativistic Landau 
levels [10| . and the Majorana equation [llj]. The goal 
of these simulations have been to experimentally observe 
analogs of long predicted but never measured effects like 
Zittcrbcwcgung and Klein's paradox. Experimental sim- 
ulation has already been demonstrated for Zitterbewe- 
gung in the ID free Dirac e quat ion 12j, and the Klein 
paradox in a linear potential [13( . 

In this Letter, we propose the simulation of neutral rel- 



ativistic particle with EDM by means of the Dirac equa- 
tion. This allows to experimentally emulate the behavior 
of EDM in an electrostatic field and two ensuing effects, 
which have not been measured so far: (i) lifting of spin 
degeneracy by an electrostatic field, and (ii) Larmor-likc 
precession of a particle spin in an electrostatic field. In 
addition we explore a few relativistic properties of these 
effects and propose their emulation with trapped ions. 

The Dirac Hamiltonian for a neutral particle with elec- 
tric dipole and magnetic moments in an external electro- 
magnetic field is given by [3l 

H = ca-p + jimc 2 + d a (0a-Bc + 2/3S-E) + 

+ Ha(0a-E/c-2pS-B) , (1) 

where c is the speed of light, m is the particle mass, a 
and $ are the Dirac matrices, S = —j6t x at is the spin 
vector operator in relativistic theory, E is the electric 
field, B is the magnetic field, d a is the electric dipole 
moment of the relativistic particle and \i a its magnetic 
dipole moment. The extended Dirac Hamiltonian (Q} is 
written in its standard form in the lab reference frame. 
For d a 7^ 0, the Hamiltonian is not invariant with re- 
spect to space inversion and time reversal 14|. For a 
neutral particle, such as neutron, neutrinos or some hy- 
pothetical dark matter particles |l| , moving in a constant 
electrostatic field with respect to the lab reference frame 
we have B = 0, which simplifies the Hamiltonian to the 
block matrix form 



H 



mc 2 l + d a cr-E cap + i/j, a a- E/c 
ca-p — in a &-E/c —mc 2 \ — d a a-E 



, (2) 



where a is a vector of the Pauli matrices and 1 is the 
2x2 identity matrix. In this fully relativistic Hamilto- 
nian the term \[i a a-Ejc coupling the magnetic dipole 
moment ^ a to the electrostatic field E in the lab frame 
of reference bears some resemblance to the spin-orbit 
coupling term which appears in the nonrelativistic limit 
of the Dirac equation. Classically the spin-orbit cou- 
pling term is induced from the Joules-Bernoulli equa- 
tions as proportional to B± ~ ox£, where v is the 
speed of the particle. However the term \[i a a ■ E j c from 
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the fully relativistic Hamiltonian is not the spin-orbit 
coupling term which although being relativistic in na- 
ture appears in the nonrelativistic Pauli equation. The 
eigenvectors of H are sought in the form of plane waves 
= \l(p))e~' p " r / n where \l(p)) are four-component 
spinors with eigenvalue Ei. The general solution of the 
time-dependent Dirac equation ifrd^/dt = is given 
by I*) = Ej bi\l(p))e-^P- r+Elt ^ h where h are complex- 
valued coefficients. 

Recently, simulation of the 3D Dirac equation in its 
supcrsymmetric representation without external poten- 
tial has been proposed for simulation of Zitterbcwegung 
of a free electron Q. Following this work we propose 
here how to simulate Eq. ^ in the standard represen- 
tation. The simulation of the 3D Dirac equation re- 
quires a single ion trapped in a Paul trap, in which the 
ion oscillates in the three spatial directions x, y, z with 
frequencies Uj (j = x,y,z). The Dirac bispinor \l(p)) 
is implemented as a linear combination of four inter- 
nal ion levels |a), |6), |c) and \d), which represent the 
internal degrees of freedom of the relativistic particle: 
\l(p)) — u a \a)+Ub\b)+u c \c)+Ud\d) . The motional degrees 
of freedom of the simulated particle can be mapped to the 
ion vibrations using the relations pj = \h{a}- — a,j)/(2Aj) 

with j = x,y,z, where a,j and Oj are the phonon cre- 
ation and annihilation operators, Aj = ^jTi/2Mvj is the 
spread in position of the ground-state wave function and 
M is the ion mass[l5j|. 

The Hamiltonian can be implemented by simul- 
taneous application of detuned red-sideband (Jaynes- 
Cummings, JC), blue-sideband (anti-Jayncs-Cummings, 
AJC) and carrier interactions between appropriately cho- 
sen pairs of the ion levels | a) , \b), | c) , | d) . The detuned 
JC and AJC Hamiltonians read 



rAJC 



a a-e 



Hf 3C = h? lj tt j (a + 



4* 



') + hSja z , (3a) 
") + h6 j a z , (3b) 



where <p r and 4>b are the red- and blue-sideband phases, 
6j is the detuning, flj are the Rabi frequencies, tr + and 
<t~ are the raising and lowering operators between two 
pairs of internal ion levels and rjj = k^Jh/2Mvj is the 
Lamb-Dicke parameter, where k is the wavenumber of the 
driving field [15|- Homogeneity of space requires to set 
the trap frequencies in the three spatial directions equal 
to each other v x = v v = 



y = v z . This ensures that f2j 
A and r/j = r/ for all spatial directions j. 



n. 



In our proposal the mass term /3moc 2 is implemented 
by the Stark shift parts of two JC and two AJC interac- 
tions applied simultaneously on the transitions \a) <R- \d) 
and \b) <-> |c), i.e., we have the mapping (3m c 2 — > 
2h8a^ d + 2Ti8a b z c . The momentum term col ■ p maps to 
the three terms 2r]AQ(af + a b x c )p x , 2r]An(af - a b y c )p v , 

2r]Art(& x c — a bd )p z , where the superscripts in the Pauli 
matrices indicate the internal ion levels between which 
the couplings arc applied. 



The term describing the interaction between the 
EDM and the electrostatic field 2d a (3S ■ E in the Dirac 
Hamiltonian ([T]) is implemented by the carrier interac- 



tion H. 



■d) 



with Rabi frequency 



Qj ' , using the mapping 2d a fiS ■ E 



2hn^>\a. 



ab 



i.cd 



er™j. The term describing the 
coupling of the electrostatic field to the MDM is im- 
plemented with the help of carrier interaction with 
Rabi frequency fr , using the mapping ifi a f3a- E/c — > 



2hn^ \-af 



•ac] 
V 



can be implemented by setting appropriately the phases 
of the carrier interactions. This establishes the following 
relationship between the parameters of the Hamiltonian 
([2]) and the trapped ion: 



(M«/c)Ej = 2HQ 
c = 2r]Ah, 



(2) 



d a Ej 



mc 



= 2TiVL 
2hS. 



(1) 



(4a) 
(4b) 



While it is possible to implement the Hamiltonian in 
supersymmetric representation this has the drawback of 
establishing correlation between the sizes of the emulated 
mass term and the emulated EDM term. The presented 
proposal for emulation of Eq. ([2]) has the advantage of 
independent experimental control of the emulated mass 
and the EDM terms. 

The essential features and effects caused by the inter- 
action of EDM with an electrostatic field remain in the 
ID limit, which is experimentally more feasible. Follow- 
ing Greiner [lif we consider the ID Dirac equation where 
the state of the relativistic particle is described by four- 
component spinors. This allows for clearer separation of 
the spin degrees of freedom from the negative and posi- 
tive energy solutions and the simulation of lifting of spin 
degeneracy and Larmor-like precession due to the inter- 
action of EDM with an electrostatic field. The use of the 
two-component form of the Dirac equation Q would not 
allow this. The ID limit of Eq. © is 

-Hid = ca x p x + f3mc 2 + 2d a /3S x E x + i(fi a /c)/3a x E x . (5) 

Note that the terms 2d a fiS x E x and i(fi a /c)f3a x E x are in- 
variant with respect to Lorentz boosts in the x-direction 
and therefore there is no effective magnetic field seen by 
the particle in its stationary reference frame. This is an- 
other point which differentiates the i(^i a /c)$a x E x term 
from the spin-orbit coupling or vx E term which will dis- 
appear in ID limit. Equation (|5|) maps to the trapped-ion 
Hamiltonian 

H 1B = 2r)Ah(af + a b x c )p x + 2U{af + a b z c ) + 

+ 2%rt 1 \af - a c x d ) - 2hfl^(af + a bc ). (6) 

The momentum terms and the mass term in Eq. ([6]) can 
be implemented by applying simultaneously two pairs 
of detuned AJC and JC interactions on the transitions 
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| a) <H- \d) and \b) O \c) with 0,. = 3ir/2 and b = it/2, 



The splitting AE 1 = E% - E\ 



-(El 



E+) is 



2rflAafp x + 2Maf = i?j c ' ad + Hi X3C ' ad , (7a) 
2i$lAo b x c p x + 2hSa b z c = H 3 x C ' bc + H^ iC ' bc , (7b) 

with ^ = 5. The implementation of the EDM term 
2d a (3S x E x requires two carrier interactions with Rabi fre- 
quency f^ 1 ) on the transitions \a) <-> |£>) and |c) |rf): 



2^( 1 )(«rf - 5^) = fl?^ = 0) + 5^(0 = tt). Us- 
ing independently and simultaneously two more carrier 
interactions with Rabi frequency fi( 2 ) on the transitions 
| a) <-> |d) and |6) <-> |c), the MDM term i(fi a /c)/3a x E x 
can be implemented: -2hn^(af + a bc ) = H' ( 



c(2). 



7r/2)+^(0 = 7r/2). 

The presence of the EDM term 2d a (3S x ■ E x in the 
Hamiltonian ([2]) causes two interlinked effects. First it 
leads to lifting of spin degeneracy in the spectrum of the 
Hamiltonian @, which is caused just by the electrostatic 
field. This effect is technically similar to the lifting of 
spin degeneracy by a static magnetic field, which under- 
pins the anomalous Zeeman and Paschen-Back effects in 
atomic physics fl7| . This is explained naturally because 
the electric dipole moment of any particle is aligned with 
its spin. However conceptually it is different from SU(2) 
breaking by magnetic field since lifting of spin degener- 
acy due to magnetic field involves only the breaking of 
time-reversal invariance, while lifting of spin-degeneracy 
by EDM - electric field coupling breaks both time and 
space reversal symmetries. An alternative explanation 
18| can be sought in the fact that the electric dipole 
term 2d a j3S x ■ E x breaks the space inversion symmetry 
[Til ] while preserving the translational invariance in an 
electrostatic field. The second effect is a consequence of 
the first: the precession of EDM and the associated spin 
around the electrostatic field, similar to the Larmor pre- 
cession of a particle spin around a static magnetic field. 

While these two effects can also be modeled within 
non-relativistic limit the simulation of the fD Dirac 
equation presents the opportunity to study two unusual 
purely relativistic features of spin-splitting by electro- 
static field. The first is the disappearance of spin- 
splitting when the mass in the Dirac equation tends to 
zero. The second is the reduction in the size of the spin- 
splitting caused by the coupling of the electrostatic field 
to the MDM embodied in the term i(/j, a /c)$a x E x . The 
mathematical analysis follows. 

For a free Dirac Hamiltonian (E x = 0) the positive 
and negative energy eigenvalues E± = ±y c 2 p 2 + m 2 c 4 
are doubly degenerate reflecting the spin degeneracy. For 
a nonzero electrostatic field (E x ^ 0) there is no degen- 
eracy in the spectrum of the Hamiltonian ([2]). 



El = ±Vc 2 Pl + E 2 x (p a /c) 2 + (mc 2 + E~daf, (8a) 
Ei = ±Vc 2 Pl + El(p a /cf + (mc 2 - E x d a ) 2 . (8b) 



AE = y/<?p2 + El^ a /cf + (mc 2 + E x d a ) 2 

-^Jc 2 p 2 x + El(pjc) 2 + (mc 2 - E x d a ) 2 . (9) 

First, in the limit m = the spin splitting vanishes, 
AE = 0, despite the fact that the EDM and electro- 
static field are nonzero. This is in stark contrast to non- 
relativistic model where the spin-splitting does not de- 
pend in any way on the mass term. Second, the term 
i(fj, a /c)$a x E x does not lead to lifting of the degeneracy 
as AE = for d a = and fx a ^ and E x ^ by it- 
self. Estimates of the EDM and MDM of several electri- 
cally neutral particles show that mc 2 3> E x d a , (fj, a /c)E x . 
Expanding Eq. © in Taylor series with respect to the 
small variable (/j, a /c)E x and making the approximations 
E\_ ~ E+ ~ mc 2 we get for the energy splitting to sec- 
ond order AE 2E x d a - ^cd a E x . Therefore 
the effect of the term i(/j, a /c)f3a x E x is to decrease the 
energy splitting caused by the EDM - E x coupling term. 
This consequence of i(p, a /c)(3a x E x again stresses its dis- 
tinction from the spin-orbit coupling term which has the 
completely different effect to lift spin-degeneracy by itself 
and thus to increase the size of spin-splitting on top of 
of this from the EDM-E X coupling. Furthermore again 
the mass of the particle m plays role in the spin-splitting 
by being one of the determinants of its splitting. The 
amount of the decrease of AE due to the MDM - E x 

coupling term depends on the ratio A = ( ^^jy^ 2 

For practically achievable strengths of E x , A < f0 -30 , 
thus for conventional experimental setups its effect can 
be neglected and the energy splitting can be written in 
the form AE = \J c 2 p 2 + (mc 2 + E x d a ) 2 
— y/ 1 c 2 p 2 x + (mc 2 — E x d a ) 2 ■ However the effect of 
i(li a /c)(3a x E x on AE can be emulated and explored in 
an ion trap setup. 

When i(p a /c)(3a x E x is neglected, the eigenspinors 
corresponding to E± and E± are given by |± \) = 
(4,-4,-l,l) T /7Vj and |± t) = ( 4,4,M ) T M, 
with Ni = ^2 + 2\vi\ 2 , Nl = ^2 + 2\vl\ 2 , v X ± = 
w i -j- y/\ _|_ W 12^ and v± = vA ± V 1 + w^ 2 , where uA = 
(mc 2 - E x d a )/(cp x ) and w t = (mc 2 + E x d a )/(cp x ). The 
expectation values of the components S y and S z are zero 
for all four eigenspinors. The expectation value of S x for 
the four-component spinors |± is +i while the one for 
|± 1) is — i confirming the interpretation of Eqs. (f5aj) 
and (|8bp as lifted spin degeneracy. 

Consider an initial state which is a linear combi- 
nation of positive energy solutions: \^f(t = 0)) = 
e -ip xX /h( b t + 1 + |) + I)). Then the general solu- 

tion of the time-dependent problem will not involve neg- 
ative energy eigenfunctions: \^(t)) = e~ iPmX / h [b i + \+ f 

)e~* E l t + b\_ |+ l)e~^ B + t ]. The expectation value of the 
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spin component Sj with respect to is 



(S^t)) = |&+| 2 (t +\Sj\+ t) + K\ 2 (i +15,1+ I) + 



I |2 / 



2Re 



(birbi(t+\s j \+i)e h 



(10) 



The overall behavior of (S) is a precession around 
with an angular frequency w, similar to the Larmor pre- 
cession of a spin in an external magnetic field. 

The experimental signature of the lifting of spin degen- 
eracy by an external electrostatic field is the Larmor-likc 
precession of the spin expectation value. Using the map- 
ping between the relativistic particle and the ion trap 
parameters the simulated precession frequency becomes 



u = 2J v 2 A 2 n 2 p 2 Jh 2 + {8 + nw) 2 



2J v 2 A 2 n 2 P 2 x /n 2 + (s- nwy 



(ii) 



The simulation requires initialization of the trapped 
ion in state = 0)) = e-' l P* x / h {b\\ + f) + b^\+ |)). 
The construction of the initial state can be done using 
the same toolbox which is employed for the simulation 
of the Hamiltonian ^ The system should first be 

cooled to its ground state. Then the motional degrees 
of freedom should be excited representing certain val- 
ues of the simulated momentum through the relationship 
Px = — o-x)/ (2A X ). Then using combination of car- 
rier interactions with appropriate timing for the desired 
values of the parameters SI, 77, A, 17, 5 one can populate 
the four ionic levels with the required probabilities. 

The dynamics of the system can be driven by the ap- 
plication of the described combination of AJC, JC stem- 
ming from a single bichromatic source and carrier inter- 
actions will simulate the Dirac Hamiltonian ©. Con- 
trolling the transitions between four energy levels is not 
trivial, but certainly feasible. Current experiments with 
trapped ions involve even more than four levels, for op- 
tical pumping, storage, ancillas, ionization etc. In the 
proposed experimental implementation we need two pairs 
of JC and AJC laser fields for the implementation of the 
momentum term, which couples levels |o) and \d), and \b) 
and |c). One can always use energy levels that form tran- 
sitions with different carrier frequencies, and/or transi- 
tions driven by different polarizations, thereby eliminat- 
ing the possibility for interference. We note that light 
shifts are not essential for the required light intensities; 
in any case they can easily be accounted for and/or com- 
pensated. These arguments arc valid also for the imple- 
mentation of the other terms in the Dirac equation. 

The dynamics of the system is manifested in a preces- 
sion of the emulated relativistic particle spin embodied in 
oscillations of the relative phase between the two eigen- 
spinors |+ f) and |+ \) with frequency w. The two eigen- 
spinors | + f) and | + \) map to the four internal levels of 
the trapped ion and the emulated dynamics of EDM pre- 
cession maps to periodic population transfer between the 



four internal ion levels with frequency u = (-ET —E\)/K. 
This frequency u> which is the signature of the emulated 
dynamics can be measured by standard ion trap tech- 
nology such as a electron shelving from any one of the 
internal ion levels. 

Supposing a realistic electrostatic field E x = 10 
MV/cm a neutron with experimentally set upper EDM 
value of around d n ~ 10~ 26 e cm will lead to spin-splitting 
AE r~j 10~ 19 eV corresponding to precession frequency 
of ui = 10~ 4 Hz; for a neutron with SM predicted EDM 
value d n ~ 10~ 32 e cm the corresponding spin-splitting 
AE ~ 10~ 25 eV and ui = 10~ 10 Hz the precession fre- 
quency is so small that it would take of the order of 300 
years for one full precession of the spin. These values put 
considerable challenge to present and future conventional 
experiments. The emulation of the Dirac equation with 
the EDM term, Eq. provides the possibility for em- 
ulation of the discussed effects, since they allow for em- 
ulated precession frequencies in the range oj ~ 10 — 10 7 
Hz. 

In conclusion, we have proposed a scheme for simu- 
lating the EDM of neutral relativistic particles within 
Dirac theory in ion traps. We have described the lifting 
of spin-degeneracy caused by an electrostatic field for a 
particle possessing EDM, and the consequent Larmor- 
like precession of the particle spin. We have predicted a 
few unusual relativistic features of the considered effects 
and proposed how they can be emulated in an ion trap. 
Furthermore this can serve as a stepping stone towards 
more involved experimental studies of the physics of com- 
bined space inversion and time-reversal violation as well 
as CP-violation. 
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and DMU03/107. 



[1] 

m 
[3; 

[4 
[5 
[6 
[r 
[8 
[9 
[10 



K. Sigurdson, M. Dor an, A. Kurylov, R. R. Caldwell 
and M. Kamionkowski, |Phys. Rev. D 70, 083501 (2004J 
M. Pospelov and ~K~. Ritz, 

Annals of Physics 318, 119 (2005)] 
E. D. Commins, S. B. Ross, D. DeMille, and B. C. Regan, 
|Phys. Rev. A 50, 2960 (1994)] 

and N. 



J. H. Smith, E. M. Pur cell 
|Phys. Rev. 108, 120 (1957) 



F. Ramsey, 



J. J. Hudson, D. M. Kara, I. J. Smallman, B. E. Sauer, 
M. R. Tarbutt, and E. A. Hinds, Nature 473, 493 (2011). 
C. A. Baker, D. D. Doyle, and et. al., 
|Phys. Rev. Lett. 97, 131801 (2006) | 
TT. Lamata, T. Leon, T\ Schatz, and E. Solano, 
|Phys. Rev. Lett. 98 253005 (2007) | 
A. Bermudez, M. A. Martin-Delgado, and E. Solano, 
|Phys. Rev. A 76, 041801 (2007) | 

J. Casanova, J. J. Garcfa-Ripoll, R. Gerritsma, C. F. 



Roos, and E. Solano, |Phys. Rev. A 82, 020101 (2010) 
A. Bermudez, M. A. Martin-Delgado, and E. Solano, 



5 



|Phys. Rev. Lett. 99, 123602 (2007)] 

[11] J. Casanova, C. Sabin, J. Leon, I. L. Egusquiza, R. Ger- 
ritsma, C. F. Roos, J. J. Garcia-Ripoll, and E. Solano, 
|Phys. Rev. X 1 021018 (2011)| 

[12] R. Gerritsma, G. Kirchmair, F. Zahrunger, E. Solano, 
R. Blatt, and C. F. Roos, Nature 463, 68 (2010). 

[13] R. Gerritsma, B. P. Lanyon, G. Kirchmair, 
F. Zahringer, C. Hempel, J. Casanova, J. J. Garcia- 
Ripoll, E. Solano, R. Blatt, and C. F. Roos, 
Phys. Rev. Lett. 106, 060503 (2011) 



[14] B. Thaller, The Dirac Equation (Springer- Verlag, 1992). 
[15] D. Leibfried, R. Blatt, C. Monr oe, and D. Wineland, 

|Rev. Mod. Phys. 75, 281 (2003) | 
[16] W. Greiner, Classical Electrodynamics, edited by D. A. 

Bromley (Springer, 1998). 
[17] A. Messiah, Quantum, Mechanics (Dover Publications, 

1999). 

[18] C. Kittel, Quantum Theory of Solids, 2nd ed. (John Wi- 
ley & Sons, 1987) p. 528. 



